home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / demos / oscil / oscil.sci < prev   
Text File  |  1999-09-16  |  649b  |  36 lines
  1. function [yd]=fty(t,y)
  2. //
  3. //  RHS of Differential equation
  4. //
  5. en = ener(t)
  6. if y(1)<0 then
  7.   if en<.5 then
  8.     yd(1) = 0;
  9.   else
  10.     yd(1) = l.*(e-r*y(2));
  11.   end;
  12.   yd(2)=-y(2)/(r*c);
  13. else
  14.   yd(1)=l*(e*en-y(2)*r);
  15.   yd(2)=(y(1)-y(2))/(r*c);
  16. end;
  17.  
  18. function [en]=ener(t)
  19. //
  20. // Energy
  21. //
  22. t1=t-pe1*int(t/pe1);
  23. if t1>pe2 then en=0; else en=1; end;
  24.  
  25. function [en]=ene1(t,pe1,pe2)
  26. t1=t-pe1*int(t/pe1)
  27. en=round( ones(t1)-t1/pe1 + (.5-(1.-pe2/pe1))*ones(t1) )
  28.  
  29. function [en]=ene2(t,pe1,pe2)
  30. t1=t/pe1-int(t/pe1)
  31. en=round( ones(t1)-t1 + (.5-(1-pe2/pe1))*ones(t1) )
  32.  
  33. function [en]=ene3(t,pe1,pe2)
  34. en=round( ones(t)-t/pe1+int(t/pe1) + (.5-(1-pe2/pe1))*ones(t) )
  35.  
  36.